Characterization of two-qubit perfect entanglers 
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Here we consider perfect entanglers from another perspective. It is shown that there are some special perfect 
entanglers which can maximally entangle a. full product basis. We have explicitly constructed a one-parameter 
family of such entanglers together with the proper product basis that they maximally entangle. This special 
family of perfect entanglers contains some well-known operators such as CNOT and DCNOT, but not y^SWAP. 
In addition, it is shown that all perfect entanglers with entangling power equal to the maximal value, | , are also 
special perfect entanglers. It is proved that the one-parameter family is the only possible set of special perfect 
entanglers. Also we provide an analytic way to implement any arbitrary two-qubit gate, given a proper special 
perfect entangler supplemented with single-qubit gates. Such these gates are shown to provide a minimum 
universal gate construction in that just two of them are necessary and sufficient in implementation of a generic 
two-qubit gate. 

PACS numbers: 03.67.-a, 03.65.Ud, 89.70,+c 



I. INTRODUCTION 

Entanglement has been proved to be a crucial ingredient in 
many quantum information processing (QIP) tasks, such as 
quantum computation and quantum communication j]]]. In 
this respect, entanglement is a unique quantum mechanical 
resource which its production, quantification, and manipula- 
tion are of paramount importance in QIR A fundamental rel- 
evant question here is how to characterize entangling capa- 
bilities of quantum operations. In fact, in this regard a lot 
of investigation have been done from many different aspects 

siiiiiaiaiiaiiiiiiaiiiiiiiiiiiaiaiii. 

In this paper we investigate the problem of characterization 
of perfect entanglers 1 17] from a new perspective. Perfect 
entanglers are defined as unitary operators that can generate 
maximally entangled states from some suitably chosen sepa- 
rable states. These operators are, in some sense, half of all 
nonlocal two-qubit unitary operators. Thus characterization 
of their (geometric and algebraic) structure(s) will be impor- 
tant. Fortunately, there are some well-defined (and relatively 
easy to calculate) tools to measure entangling properties of 
quantum operators, and also entanglement properties of two- 
qubit states 12 0, fl^l . By using these tools we show that 
one can analyze perfect entanglers and find that some of them 
have the unique property of maximally entangling a complete 
set of orthonormal product vectors. We explicitly construct 
such operators and show that they have the maximal entan- 
gling power, I . It is proved that these are the only unitary op- 
erators which have this property. However this investigation 
has its own importance, we provide a probable application for 
them. An analytic way to use such special perfect entanglers 
as the building blocks of the most efficient universal gate sim- 
ulation is put forward. This assigns another remarkable role 
for such special entanglers in universality issues and the re- 
lated topics. 

The structure of this paper is as follows. In Sec. |H] we 



review some necessary definitions and some important facts 
about nonlocal two-qubit operators. This section is supple- 
mented with appendices at the end the paper. In Sec. 
we explicitly construct a one-parameter family of such spe- 
cial perfect entanglers and some relevant examples are pro- 
vided. Also we go further and show that the family we have 
found is exactly the family of perfect entanglers with maxi- 
mal entangling power, |, and no other perfect entangler ex- 
cept these ones have the mentioned property. In Sec. lIVI we 
discuss applications of these special perfect entanglers in uni- 
versal gate simulation. We show that they can be used as el- 
ementary building blocks in minimum universal construction 
of two-qubit gates. The paper is concluded in Sec.lVl 



II. PRELIMINARIES 

In this section we want to review briefly some facts on non- 
local two-qubit unitary operators. More discussion on this 
matter can be found in the appendices and references therein. 
It is known that fl frTlEoll any two-qubit unitary operator U 6 
SU(4) can be written in the following form 

U = (A 1 <E) B 1 )e-' l( ' Cl ' Tl ®' Tl+C2 ' T2 ®' T2+C3 ' T3 ®' 73 ')(A2 ® B 2 ),(\) 

where oVs are the Pauli spin matrices, and At, Bi e SU(2) 
are some single-qubit unitary operators. By using this decom- 
position for U, called the canonical decomposition, one can 
define local equivalence of two unitary operators. Two opera- 
tors U, V E SU(4) are considered locally equivalent, U = U', 
if there exist some single-qubit unitary operators ui, U2, «i, v 2 
€ SU(2), such that 



U 1 = ui ® v\Uu 2 <8> v 2 . 



(2) 



•Electronic address: tayefehr@mehr.sharif.edu 



It has been shown that using the concept of local equivalence 
one can find some freedoms in Cj's such that one can always 
restrict oneself to the region f > c\ > c 2 > | C3 1 in the 
(ci, c 2 , C3) space (Appendix |AJ. This is the so-called Weyl 
chamber 11711 . Hereafter, whenever we write [01,02,03] we 
will assume using these freedoms. So when we want to deal 
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with local equivalence we can always consider U = U = 
[ci, c 2 , C3], up to the freedoms. For the following uses it is 
helpful to present the explicit form of U. In the standard com- 
putational basis (where (T3IO) = |0), cr 3 1 1) = —|1)), we have 

/ e~ lC3 c- -ie- lC3 s- \ 

-ie lC3 s + e lC3 c+ A ' 

\ -ie~ lC3 s- e- lC3 c~ j 

where c = cos(ci ± c 2 ) and s = sin(ci ± c 2 ). Consider 
the usual Bell basis, {|$±) = JMglll , |*±) = 
Then, it is easily seen that in the magic basis 1 

|*i) = |*+), |* 2 ) = -i|4->, (4) 
|3 3 > = |*->, |$ 4 ) = -i|¥+>, (5) 

W is diagonal 

W = Et= 1 e- <A *|*fcX**l. (6) 

where 

Ai = ci - c 2 + C3, A 2 = — ci + c 2 + C3 
A3 = — (ci + c 2 + c 3 ), A 4 = ci + c 2 - c 3 . 

Makhlin 1 9] and Zhang et al. [ IT] have shown that the local 
equivalence classes can be characterized uniquely by the two 
local invariants G\\U] and G 2 [C/] which are as follows 

Gi = i[e- 24C3 cos2( Cl -c 2 ) + e 2jC3 cos2( Cl + c 2 )] 2 , (7) 
G 2 = cos(4cx) + cos(4c 2 ) + cos(4c 3 ). (8) 

Moreover, in l9 Hl7ll the important notion of a perfect entan- 
gler has been defined. Strictly speaking, a two-qubit unitary 
operator is called a perfect entangler if it can produce a maxi- 
mally entangled state from a suitable unentangled (separable) 
one. It is not hard to see that [CNOT], [ DCNOT] 2 , and [Vswap] 
classes are all perfect entanglers. For example, for [CNOT] a 
choice can be [CNOT]|0)|0) = |00> ^l |n> . 

To characterize entangling capabilities of unitary operators 
entangling power has been introduced 1 2 ] . This measure is de- 
fined to be the average entanglement that the unitary operator 
can produce when acting on separable states. It can be shown 
that entangling power of a unitary operator U, e p (U), has the 
following simple form 

e p(U) = e p (li) = j%[3 — (cos4ci cos4c 2 

+ cos 4c 2 cos 4c3 + cos 4c3 cos 4ci )] , (9) 



1 There are different conventions about the magic basis I5 ll3tll9ll . We fol- 
low the definition in 

2 DCNOT or double-CNOT is defined as DCNOT=CNOT 2 iCNOTi2, 
where the subscripts 12 for CNOT mean that the first (second) qubit is the 
control (target) qubit, and similarly for 21. 



which only depends on the nonlocal part [01,02,03] (Ap- 
nendixlBl. 

Other important points that are helpful in our next discus- 
sions are separability and maximally entanglement conditions 
in the magic basis. It has been shown in |3[] that a state 
I*) = Efe Mfc|$fc) is separable iff J2k /"J = °- As well, |#) is 
maximally entangled iff /1 2 = e~ lS \^k\' 2 (k = 1, . . . , 4). This 
means that a state is maximally entangled iff its coefficients in 
the magic basis are real, up to a global phase. 

III. SPECIAL PERFECT ENTANGLERS 

In this section, we are going to study some special members 
of the space of perfect entanglers. As we see below this turns 
out to be useful in characterization of structure of the space 
of perfect entanglers. In this regard, we pose the following 
question: 

Does there exist any perfect entangler which can maximally 
entangle a full separable basis? 

This question is motivated by the following observation which 
was firstly reported in 112 In . To be complete in our discussion 
we reproduce it here. It can be seen that for the [CNOT] class 
one can find four orthonormal product states which are trans- 
formed to some maximally entangled states. For example, we 
have 

[CNOT]|00) = ^ 7 p^, [CNOT]|01) = MM 

[CNOT]|10) = M^i>, [cNOT]|ll) = Jli«). (10) 

This property may not seem so important at this stage. But if 
we want to find such a separable basis for the [^SWAP] class 
the result, as is shown below, is a failure. To see this, we 
note that the most general separable basis (up to general phase 
factors for each vector) that can be considered is as follows 

|*i) = (o|0) + 6|l))®(c|0)+d|l)), 

|*a) = (-6|0) + o|l))®(c|0) + d|l)) ) 

|* 3 > = (e|0)+/|l))®(-dl0)+c|l)), 

|* 4 > = (-/]0)+e|l))®(-d]0)+c|l>), (11) 

in which bar means complex conjugation and |a| 2 + \b\ 2 = 
|c| 2 + \d\ 2 = |e| 2 + |/| 2 = 1. Now if we act on one of these 
states, say |\f'i), by the [SWAP] class we obtain 

[yswAP]|*i) = e- l i ac|00) + ^{ad - i&c)|01) 

-iZjL(ad + ibc)\10) + e- i %bd\ll). (12) 

This state is maximally entangled if its concurrence 1 19j, C = 
\ad-bc\ 2 , is equal to 1 (Appendix|SJ. Thus the [Vswap] class 
turns these orthonormal states into maximally entangled ones 
if the following conditions are fulfilled simultaneously 

\ad-bc\ 2 = l, \ac+bd\ 2 = 1, (13) 
|c/-de| 2 = l, |ce + J/| 2 = l. (14) 

However, by side by side addition of the equations in 1131 . 
we get (|a| 2 + |6| 2 )(|c| 2 + |d| 2 ) = 2 which in the light of 
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the normalization conditions is an apparent contradiction! A 
similar contradiction is obtained from Eq. (I14> . Therefore the 
[v / swap] class cannot maximally entangle any full product 
basis. However, it must be noted that this fact does not forbid 
finding a pair of orthonormal states which become maximally 
entangled by this class. For example, by choosing the follow- 
ing pair of orthonormal product vectors 

= (o|0) + 6|l))®(c|0) + d|l)), 
I*;) = (-6|0) + fi|l» <g> (-J|0) + c|l)), (15) 

in which \ad—bc\ 2 = 1, no contradiction is turned out. Itis in- 
teresting to note that e p (^/sWAP) = | which is the minimum 
entangling power among all perfect entanglers. 

In the following we are going to find some other unitary 
operators having this property. Hereafter, the perfect entangler 
that can maximally entangle a full product basis are called 
special perfect entanglers (SPE). We are trying to explicitly 
construct such operators together with their related product 
basis. 

Since working with a general product basis, Eq. il It . is hard 
let us, for the time being, restrict ourselves to a rather special 
case of product basis set. It is seen that by a suitable local 
unitary operator, U\ ® U2, it is always possible to transform 
two of the basis vectors of {|^i)}f = i, say \^>i) and \^ 2 ), to 
the separable states |00) and |10). If we choose 

"-(iH).*-^?)- <16) 

then U\ ® U 2 acts on the vectors of Eq. il It as 

C7i®tf 3 |*i) = |00), 
U 1 ®U 2 \t>2) = |io), 

Ui ® I7 2 |*3> = [(fie + 6/)|0) + (-be + a/)|l)] 
UxZsU^i) = [-(-be + af)\0) + (bf + ae)\l)] |1). 

(17) 



The maximal entanglement condition for these states requires 
that 



\e 2iC3 c+8+A 2 - e- 2iC3 c~s-B 2 \ = 
\ e 2ic 3c - s - A 2 _ e -2ic 3c + s + S 2| = 1 (19) 

As is seen the conditions on c\ and c 2 are independent of the 
the last conditions. Combination of this set of equations gives 
the following results 

sin2(ci -02) =±1, sin2(ci + c 2 ) = ±1, (20) 
\e 2iC3 (±l)A 2 - e- 2lC3 (±l)S 2 | = 1. (21) 

By using the triangular inequality, \zi + z 2 \ < \zi\ + \z2\, for 
two complex variables z% = |^j|e 8 * >i , and noting that the equal- 
ity holds when <pi — <p 2 (mod 2tt), Eq. ( I21> can be simplified 
more. If we note that \e 2iC3 A 2 ±e" 2iC3 B 2 \ < \A\ 2 + \B\ 2 = 1, 
and parameterize A — cos 8 and B = e 2 ^ sin 6 (2<p, as the 
phase, is chosen for later convenience), then \e 2lC3 cos 2 9 ± 
e -2i( C3 -20) sm ^| _ 1 requires that 

g2'ic 3 _ _|_ e -2i(c 3 -20) 

which has the following simple solutions 

Now if we use the freedoms in choosing Cj's for a class (Ap- 
pendix|X|i, there only remain the following solutions 

c 3 = 0, f + <P. (22) 

We should now go back to the condition J20l i. Apparently four 
different cases can be obtained from this equation. But again 
by a simple algebra and using the freedoms in Cj's of a class it 
appears that the two classes 



Thus this separable basis is locally equivalent to the following 
more simpler set 

{|00>, |10>, (A|0) + S|l»|l>, (-B|0) + (18) 

for the specific values of A and B as in Eq. J17I . Here a point 
is in order. However for any unitary operators Ui,U 2 € SU(2) 
we have E{Ui®U 2 U\$>)) = E(U\&)), generally it cannot be 
concluded that E(JJUi®U2\^)) = E(U\^)). As an example, 

cnot|00> = |00), but cnot(H ® 1)|00) = |00) ^1 11) (H is 

the Hadamard matrix). This point hints us that our solution(s) 
may not be the most general one(s), that is, it is still possible to 
find other SPE in addition to what we find below. We return 
to this point a bit later. Anyway, action of a U on the basis 
vectors in Eq. dl8> . respectively, gives the following results 



\ 






-ie~ iC3 s- 
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■B\ 



I 



e iC3 c+A 
-ie iC3 s+A 
e- iC3 c~B ] 



-ie~ lC3 s'A \ 
-e iC3 c+B 
ie iC3 s+B 
e- iC3 c-A j 



Ui 



u 2 



[o,f,< 

1 - 



(23) 
(24) 



are our independent solutions. Furthermore, by a simple re- 
definition of <p in Eq. J24> and another use of the freedom 
for Cj's, it can be inferred that we actually have found a one- 
parameter family of SPE 



,0], O<0< 



(25) 



Figure Q shows the geometrical representation of the SPE 
family in the (01,02,03) space 3 . It should be noted that in 
the two-qubit case we are fortunate to have this simple geo- 
metrical interpretation, which is relatively absent for higher 
cases. 



' The SPE line A1A2 in Fig.Qcorresponds to the A2L line in Fig. 1 of |34 
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Now, taking into account the last points of Sec. [II], separability 

of {l^sep)}^! and maximally entanglement of {l^me)}^!) 
respectively, requires that 



E fe /4 i} = o, 

from which (by using Eq. &291 ) one obtains 



FIG. 1: Geometrical representation of SPE gates on the (ci, Ca, C3) 
space. The line A\Ai is the location of all possible SPE. The points 
Ai, B, and A 2 denote the [CNOT], the [B], and [DCNOT] gates, re- 
spectively. The filled area is projection of the Weyl chamber which 
is the region containing nonlocal gates. 



Using the form, for example, in Eq. d23l we can simpl^ 
find the form of the transformed states (up to local unitary 
operators), which are as 

--i^ |00)-H|U) r ij, -i|01) + |10) 

^[sin(9(i|00) + |11)) + cos0(|Ol) - i|10»], 
^[cos0(z|OO) + |ll))+sin0(t|lO) - |01))], (26) 
and are clearly maximally entangled. 

It is interesting to see that the [CNOT] gate (for <j> = 0) and 
the [DCNOT] gate (for <f> — j) are both members of this one- 
parameter family of SPE, however, [VswapJ is not. Also the 
local invariants of the SPE in Eq. (I25> are 



0. 



(i) 1 ziitfM 1 (i) 1 

w) — ±e 2 \Wf, I 



(31) 
(32) 

(33) 
(34) 



Gi=0, G 2 =cos4</>. 



(27) 



Up to now, we just have shown that at least there exists 
a one-parameter family of unitary operators that are SPE, 
Eq. d25i . Now, we consider the problem in a general man- 
ner. In fact, we show that the family J25I is the only possible 
family of operators that are SPE. To this end, we use the di- 
agonal form of U in the magic basis, Eq. (|6}. In this basis we 
can write 



I*. 



where {|* s ( $)}f =1 is a general product (separable) basis, and 

similarly {|^me)}i_x * s a maximally entangled basis. So our 
question, now, is to find unitary operators, U, for which we 

(i) 

can find some suitable separable basis (fi k 's) and a maxi- 
mally entangled basis (w^) such that 

«l*g> = l*S2>- (28) 
By using Eq. the following relation between fij^ 's and 



's is obtained 



w 



(0 



-iA fc ,,W 



a4 (*,* = !,.. .,4). 



(29) 



In addition, orthonormality conditions, (^senl^se ^ — 
) = «« read as 



- sep I * sep / 



Erf/^E^' 



(30) 



If we replace Eq. (I34t in Eq. J33l > and multiply it by e 2iC3 , after 
separating real and imaginary parts we get 

| 2 + \wP | 2 cos a 2 + \w { 2 i} | 2 cos a 3 + \w[ i} \ 2 cos a Y = 0, (35) 
\w[ l) \ 2 sina 2 + |w 2 ^| 2 sin 0-3 + \w^\ 2 sinai = 0, (36) 

for i = 1, . . . , 4, where 

ai = 4(ci + c 2 ), a 2 = 4(ci + c 3 ), a 3 = 4(c 2 + c 3 ). 

Orthonormality condition, Eq. d30l . is equal to unitarity of the 
matrix Uik ■= w^, which instead gives the below results 



E fc k l) l 2 = Erfl 2 = i- 



Efe^fe 



W„„00 



( fe ) 



(fc) O) X 



(37) 
(38) 



By using Eq. (I34> . Eq. d38l > can be simplified as 

e4(' M -' W ) Efc l^rfl = E fe k^ll-fl = ^-(39) 



If we now add the relations in Eq. d35l for different i's, and 
consider Eq. (I37> . it is obtained that 



cos ai + cos a 2 + cos a 3 = — 1 . 
In a similar manner, Eq. d36l results into 

sinai + sina 2 + sina 3 = 0. 



(40) 



(41) 



By combination of these two equations, the following final 
condition for a class [a, c 2 , c 3 ] being a SPE is obtained 

cos 4ci cos 4c 2 + cos 4c 2 cos 4c 3 + cos 4c 3 cos 4ci = — 1 .(42) 

This, in turn, implies that e p (U) = |. Therefore, we have 
proved that if a unitary operator, U, is a SPE it must have the 
maximal entangling power. Now, let us investigate all possible 
solutions of Eq. d42i . It is trivial to check that if one of the Cj's, 
say ci, is equal to ?(0), then the only choices for the other c/s 
are 0( f) a nd < < j. This is the very family of SPE as 
in Eq. (125 > . More generally, the only acceptable solutions of 
Eq. (1421 are also of the form in d25i (and its permutations). 
Thus it has been shown that only the unitary operators [j , <j>, 0] 
which are the unique family with e p = | can be SPE. This fact 
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stresses that for the SPE entangling power is a sufficient tool 
which can characterize them uniquely. 

Now let us consider the subspace of non-perfect entanglers, 
especially the [SWAP] class in the sense that how they act on 
separable bases. We want to point out a simple property of 
this class. As has been proved in l22ll any nonlocal two-qubit 
(qudit) unitary operator other than those in the [SWAP] class 
can be used as a universal operator (of course together with el- 
ementary single-qubit unitary operators). Additionally, it has 
been known that any unitary operator V which leaves product 
states product (that is, for any \x) and \y) there exist some \u) 
and \v) such that V^|cc) = is either a local unitary 

operator A® B or (A® B)S, for some A and B. Now, it can 
be deduced that the only unitary operators that can transform 
any separable basis to some separable basis are those in the 
[SWAP] class. This can be seen as follows. Consider two arbi- 
trary vectors \a) and \b). Suppose that V is such an operator 
that transforms any product basis to another one. To see that 
if V|a)|&) is a separable state, we can construct another basis 
set in which \a)\b) is one of the basis vectors. For example, 



{\a)\b), \a x )\b), |c)|& J 



is such a new basis. Thus V must transform \a)\b) to some 
|a')|6'). As mentioned above, this indicates that V is ei- 
ther local or locally equivalent to the SWAP operator. Here 
we should remark that the above statement is true only if 
V leaves any separable basis a separable one. For example, 
the CNOT operator transforms the usual computational basis 
{|00), 1 01) , 1 10), |H)} to another separable basis, however, it 
is nonlocal. 

Finally we show that by using the spectral decomposition 
of the nonlocal part, Eq. 0, a simple (yet remarkable) prop- 
erty of U, which is related to the magic basis, can be deduced. 

Lemma. There is no separable basis {\oti) <£> \f3i)}t = i that 
can be transformed to the magic basis by an operator U, as in 
Eq. 0. 

Proof. (By contradiction) Suppose that for a given operator 
U there exists such a basis, U\a.i)®\f3i) = This is equiv- 
alent to W\$>i) = \cti) <g> |/3j). By using the form in Eq. (J6J, 
this gives rise to e lXi \^i) — |aj) (8 |/3j) which is an evident 
contradiction! □ 

It should be emphasized that, for example, for the [CNOT] 
class we have 



|Q)±|i) 

V2 



|0) 



|00)±|11) 
V2 



|0)±|1) l-n , |oi)±|io) 

The right hand side is the usual Bell basis. However, this does 
not violate the above lemma. In fact, the point is the difference 
between the usual Bell basis and the magic basis. The two 
sets of basis vectors are not locally related. This can be seen 
by examining, for example, the following operator 



|* 2 )($- 



|$ 4 ><*+| 



This operator is diagonal in the magic basis, and thus is non- 
local. Actually, for this specific operator Gi[J7] = and 
G2[U] = -1 which indicate U £ [DCNOT]. 



IV. APPLICATIONS TO EFFICIENT TWO-QUBIT GATE 
SIMULATION 



In this section, we are going to present a remarkable impor- 
tance of the SPE, which can put forward special theoretical 
and experimental interests in them. 

From theoretical and experimental viewpoints, efficient 
quantum gate synthesis is of high interest. In fact, there is 
a vast literature on this topic with many remarkable findings 

EH HIH M \m\M M El IB HUMS, it has 

been known that a general (multi-qubit) quantum gate can be 
simulated using a quantum circuit built of elementary gates 
which operate on single and two qubits 12311 . For instance, 
Barenco et al. |23] have shown that a combination of the 
CNOT and single-qubit gates is universal, in the sense that any 
unitary operation can be simulated by them. This is the com- 
monly adopted universal set in paradigm of quantum compu- 
tation. However, there has been a lot of interest in universal 
sets which are highly (or the most) efficient. Actually, there 
is a practical reason behind this demand. In present day ex- 
periments, two-qubit gates as the CNOT gate can only be im- 
plemented imperfectly due to today technological limitations 
(the so-called decoherence). Therefore, in order to reduce the 
probability that an error occurs in performing a certain uni- 
tary operation on several qubits, instrumentally it is highly 
demanded that the number of the building block gates (i.e., 
the number of times that the qubits interact) be as small as 
possible. Furthermore, it is clear that these attempts can be of 
advantage in analyzing the algorithmic complexity of a given 
quantum computation and other QIP related tasks 1 1 8 ] . 

Very recently, it has been shown that in order to implement 
an arbitrary two-qubit unitary operator not more than three 
applications of the CNOT or the DCNOT gates are necessary 
yUlH]. In this sense, the cnot and the dcnot gates are 
highly efficient elementary gates. However, remarkably, it has 
also recently been shown that the B = °i®<^) 
gate is the most efficient one, in that any arbitrary two-qubit 
unitary operator can be constructed only by two applications 
of such gate 1 34]. The fact we would like to stress on is that all 
of these three gates; the CNOT, the DCNOT, and the B gates, 
are SPE. This can raise the question that if other SPE gates 
also can give rise to a minimum universal gate construction. 
To this end, we follow the method proposed in mm. we 
need to remind that the B gate as the most efficient known uni- 
versal gate can simulate a generic nonlocal two-qubit operator 
with the following circuit 











B 




U 













(43) 



where the parameters a and b satisfy 



sin 2a 



COS 2C2 COS 2C3 

1 — 2 sin 2 C2 cos 2 C3 : 
2 „_ „„„2 . 



cos 26 = 1 — 4 sin c-i cos C3. 



(44) 



Equivalence of these two circuits can be proved simply by 
showing that the nonlocal invariants (G\ and G2) of the both 
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circuits are the same. For our study, we choose an arbitrary 
SPE, which we call C[0] = e - l ^ ai ®' yi+<t " y2 ® (y2 \ We are 
looking for some SPE gates other than the B = ] gate 
that give rise to a minimum construction of two-qubit unitary 
operators universally. Following 13411 . we consider the circuit 
below which is similar to the one in Eq. (I43> 



C[0] 




e-icio-2 




cm 

















(45) 



and try to find the parameters a and 6 such that this circuit can 
simulate a generic nonlocal gate U(c\, C2, C3), as in Eq. Q. 
After computing the local invariants of this circuit and com- 
paring them with the invariants of [c\, C2, C3], it is seen that 
the following equations must be satisfied 



cos 2a sin 26 sin 40 = sin 2c2 sin 2c3 (46) 
(1 + cos 40) cos 4a + (1 + cos 40) cos 4a cos 2b + (3 cos 40- 1) cos 26 + 4cos2c 2 cos2c 3 - (1 + cos 40) = (47) 

I 



In fact, there is another (relatively long) equation which also 
must be satisfied as the consistency check for the possible so- 
lutions, which we have not written here. By the way, it is easy 
to see that the following two solutions can be found 



cos 26=1 



( 1 -cos 2c 3 ) ( 1+cos 2c 2 ) 
1 — cos 40 



COS 2 2a = ( 1 + cos2c 3)( 1 - cos 2c 2 )tan 2 20 

1+cos 2C3 —cos 2c2 — COS 2C2 cos 2C3 : 



and 



cos 26=1 



(1+cos 2e 3 )(l-cos2e 2 ) 
1 — cos 40 



COS 



2 2 (1— cos 2c3)(l+cos 2c 2 ) tan 20 



(48) 



(49) 



1 —cos 2c3 +cos 2c 2 +cos 2c 2 cos 2C3 ' 

Both of these solutions satisfy the mentioned consistency 
check relation, for any 0. However, from these solutions it 
immediately can be deduced that we must have 7^ 0, f . 
This means that the CNOT and the DCNOT gates cannot be 
among these gates. This fact is in accordance to the results of 
H31L 13311 which indicate that except some special cases at least 
three applications of the CNOT and the DCNOT gates are nec- 
essary and sufficient for generic gate simulations. In addition, 
after checking that the solutions must satisfy the the relations 
— 1 < cos 26 < 1 and < cos 2 2a < 1, it is found that if we 
want to have at least one solution we must choose the value of 
by considering the values of C2 and C3 of li{c\, C2, c 3 ) op- 
erator which is to be simulated. That is, not any (and hence 
C[0]) is proper for simulation of an arbitrary nonlocal gate. 
It can directly be checked that the only value of which is 
applicable in simulation of all two-qubit unitary operators is 
= Thus, in this way, we argue that the most efficient 
gate that can be used in universal two-qubit gate construction 
is the B = C[j] gate. This fact gives a unique and special 
role to the B gate among other SPE 13411 . 

To complete our discussion, in the following we show how 
one can simulate some well-known gates, such as the CNOT 
and the DCNOT, using our method. For the case of the CNOT 
gate, we have the following equivalent circuit 



C[0] 








C[0] 














CNOT 













which is valid for < < ?. Similarly, for the DCNOT gate, 
the following circuit gives an equivalent construction 



C[0] 




e -<f "1 




C[0] 












e -*i*3 e -£!><r2 e -*i<'3 





where 



cos 26 = - cot 2 20, §<0<f. 



As the final point, we want to emphasize that how the C[0] 
gates can be simulated in the standard gate construction, i.e., 
by applications of the CNOT (or DCNOT) gate. In |3ll 13311 . 
it has been shown that any nonlocal gate as [ci, C2, 0], where 
f > ci > c 2 > and ci > (except for [f ,0,0]), can be 
simulated with only two CNOT or DCNOT gates. This means 
that all SPE, the C[0] gates (0 < < f ), can also be con- 
structed very efficiently using the CNOT or the DCNOT gates. 

In summary, what we have shown is that for optimal con- 
struction of any arbitrary two-qubit gate we can always find (at 
least) one circuit synthesis by some suitable SPE. In our an- 
alytic construction, the choice of the convenient SPE is gov- 
erned by the values of c 2 and c 3 of the gate which is to be 
simulated. From these proper SPE one can pick the one that 
is available or can be implemented easily in experiment. In 
this way, typically there are more options at hand that one can 
utilize in experimental implementations. In addition, we have 
clarified the unique property of the B gate, that is it poses the 
most efficient elementary gate suitable in all generic two-qubit 
gate constructions. 



V. CONCLUSIONS 

It has been shown that some of perfect entanglers have the 
unique property that they can maximally entangle a suitable 
full separable basis. Some well-known two-qubit gates such 
as the CNOT and the DCNOT are among such special opera- 
tors, however -\/swap is not. In our case of study entangling 
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power has appeared as a useful tool in characterization of spe- 
cial perfect entanglers. The perfect entanglers with maximum 
value of entangling power, |, have turned out to be the only 
possible special perfect entanglers. 

As a possible application of these special gates their impor- 
tance in universal two-qubit gate construction has been em- 
phasized. We have provided a more general (analytic) result 
for optimality. In this regard, it has been proved that if we have 
the SPE gates at our disposal we can simulate any generic two- 
qubit gate by only two applications of the same gates. In fact, 
in this manner we have a vast possibility for universal gate 
construction. Although not any SPE gate may be applicable 
in simulation of a specific two-qubit gate, it is always possible 
to find some proper SPE for this mean. Specially, the unique 
role of the B gate, namely, its usage in all two-qubit gate sim- 
ulations, has been pointed out. In this sense, these results can 
shed new light on theoretical and experimental investigations 
to characterize the rich structure of the space of nonlocal uni- 
tary operators and using them as the most efficient universal 
elementary gates. 
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TABLE I: Nonlocal properties of some equivalence classes. 



Operator 


_ci, C2, ca] 


Gi 


G 2 


e p ° 


CNOT 


[f ,0,0] 





1 


2 
TT 


DCNOT 


[f.f.O] 





-1 


2 
9 


B 


[f,f,0] 








2 

a 


SWAP 


4", —] 


-1 


-3 





VSWAP 


[■§ ' ¥ ' 8'] 


4 





i 

6 


controlled-t7 


[x, 0,0]* 


Gi 


2Gi + 1 - 


-cos 4k 
9 


A <g> B 


[0,0,0] 


1 


3 





"This is entangling power which is defined in Eq. IC1I. 
*The exact value of x depends on explicit form of U. For example, 
/ case \ thenx= f.andGi = cos 2 0. 

^ -,i„8 e '+ cas6 J * 


if U = 



Then define Ub '■= Q'UQ, and accordingly m(U) := U b Ub 
(T means transposition). Gi[t7] and G?2[f] defined as 



Gi[U] 
G 2 [U] 



Xi 2 [m(U)} 
16 det U ' 

Xi 2 [m(U)} -Xi[m 2 (U)} 
4det£7 : 



(B2) 
(B3) 



have been proved to be locally invariant |0, fl7il . Therefore 
they can characterize equivalence classes. In this respect we 
can write Gj[J7] = Gj[ci, C2, C3]. For a nonlocal operator, as 
in Eq. Q, a simple algebra reveals that 



APPENDIX A: POSSIBLE FREEDOMS IN NONLOCAL 
PARTS 

It can be simply checked that in a triplet [ci, C2, C3], which 
represents nonlocal part of an operator, one can easily add any 
integer multiple of ^ to every Cj . This can be verified by mul- 
tiplication of the local operators ±iai®<Ji = e ±t 2< Ti ®' yi , Also 
[c%, C2, C3] = [cj, Cj, Cfe], where (ijk) is an arbitrary permuta- 
tion of (123). This can be verified as follows. If in the identity 
e u 1 ®v 1 Vu 1 ®v 1 _ Ml (g, Vie v u ^ <gj v ^ where all operators are 
unitary, we take u\®v\ — e~' l T ai (g)e~ 4 T <Tl it is obtained that 
[ci, C2, C3] = [ci, C3, C2]. A similar procedure can be used to 
prove the general statement. The other fact about [c\ , c 2 , C3] 
is that one can always make at least two of tVs positive (or 
negative). This can be verified by taking u% «i = H<Jk <E> 1 
in the above mentioned identity. To see more on the freedoms 
refer, for example, to 1 1 311 . 



APPENDIX B: THE LOCAL INVARIANTS 



Consider the matrix Q which is 



/ 1 i 

i 1 

Q i -1 

V 1 -i 



(Bl) 



[e" 2iC3 cos2(ci - 02) + e 2lC3 cos2( Cl + c 2 )] 2 , (B4) 



G 2 = cos(4ci) + cos(4c 2 ) + cos(4c3). 



(B5) 



Table|I]summarizes the results for some examples. As is seen, 
the [CNOT], [DCNOT], and [ V / SWAP] all belong to different 
equivalence classes. There is a powerful theorem which ex- 
actly determines which unitary operators are perfect entangler 

Bo. 

Theorem: A two-qubit unitary operator is a perfect entan- 
gler iff the convex hull of the eigenvalues of m(U) contains 
zero. 

An immediate result of this theorem is that among controlled- 
U operators only the [CNOT] class is perfect entangler (Ta- 
ble|D>. 



APPENDIX C: ENTANGLING POWER 

Entangling power, e p (U), is a widely accepted measure for 
quantifying entangling capabilities of a unitary quantum op- 
erator U Ql- Consider a bipartite quantum system with state 
space of TLa®'Hb ( m qubit case Ha = 7~(.b = C 2 ). If E is an 
entanglement measure over Ha <8 Hb, the entangling power 
of U (with respect to E) is defined as 



5 P (U) := E(U\il>) ® m 



(Cl) 



where the bar denotes the average over all the product states 
\ip) (£) \<j>) distributed according to some probability density 
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p{i/), 4>) over the manifold of product states. To ensure lo- 
cal invariance for e p (U) the measure of integration, p, must 
be taken locally invariant. A simple and proper distribution 
which satisfies this requirement is the uniform distribution. In 
121 E has been taken to be linear entropy, which is defined as 



Em 



AB 



1 



tr(p 2 



MB) 



(C2) 



where Pa(b) — ^ b(A){\^) ab{^\) is the reduced density ma- 
trix of the system A(B). Now it is in order to mention the 
relation of E and the other (rather more well-known) mea- 
sure of entanglement, concurrence 1 19]. Concurrence of any 
two-qubit density matrix, pab, is defined as 



C(pab) ■= max{0, A x - A 2 - A 3 - A 4 }, 



(C3) 



in which Ai > A2 > A3 > A4 are square roots of eigenvalues 
of the matrix pabGi <8> viPabVi ® a i (* means complex con- 
jugation in the standard computational basis). For a pure state 
oftwoqubits, |#) = o|00)+6|01)+c|10>+d|ll>, this reduces 
to C(W) = 2\ad — bc\. In this case, one can see that eigenval- 

. , ■ x i±^/i-c 2 (*) 
ues of the reduced density matrices are A± = — - — 5 > 

from which 



±C 2 (*). 



(C4) 



Since C ranges from (for separable states) to 1 (for maxi- 
mally entangled states), we have < E < \. 

Some important properties of e p (U) are listed below, 
(i) For every U & SU(4) we have < e p (U) < §. 



(ii) For every A,B € SU(2), e p (A ® B) 
from the fact that EQip) ® \<f>}) = 0. 



0. This is clear 



(iii) e p (U) is locally invariant, that is, for every A, B S 
SU(2), e p (A ® BU) = e p (UA <g> B) = e p (U). This property 
is also simple to verify by noting that 



E{UA\ip) ®B\<f))) 



\i>),\<t>) 



(iv) For every U € SU(4), e p {W) = e p (U). The equality 
holds because of = J5(|*». 

(v) If S 1 is the usual SWAP operator (S\ip) <8> |</>) = |0)®|^>), 
for every and |</>)), then e p (SU) = e p (US) = e p (U). 
These are immediate consequences of E(S\*f?)) = E(\^f)), 
for all I*). 

(vi) Another trivial property is that e p (S) = 0. 

To see more about the entan gling power, its properties, and 
generalizations refer to II2L fl4l il5ll . 

By using the above facts, one can find an analytical form for 
entangling power of a unitary operator as in Eq. l|9}- The last 
equality follows from Eq. iC4i and using the parametriza- 



1 I cos |e l 2 

tion \ip) = 1 1 



(0 < 9 < tt,0 < < 2tt) 



for a single-qubit state vector. Then a simple integration as 
e p (U) = J 0} sin OdOdcf) sin 6' dO'dcf/ gives the result. 
This form for e p is clearly locally invariant, and shows that 
all unitary operators in the same equivalence class possess an 
equal entangling power, however the converse is not true gen- 
erally. It is simple to check that this expression respects all 
the freedoms mentioned earlier for Cj's of a class. The last 
column in TableUgives e p for some well-known equivalence 
classes. It is important to note that the single parameter e p 
cannot replace the two local invariants G\ and G2 to uniquely 
characterize equivalence classes. 
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